Role of quasiparticles in the growth of a trapped Bose-Einstein condensate 
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A major extension of the model of condensate growth introduced by us is made to take account of 
the evolution of the occupations of lower trap levels (quasiparticles) by scattering processes, and of 
the full Bose-Einstein formula for the occupations of higher trap levels, which are assumed to have a 
time independent occupation. The principal effect is a speedup of the growth rate by somewhat less 
than an order of magnitude, the precise value depending on the the assumptions made on scattering 
and transition rates for the quasiparticle levels. 



In a previous paper [fjj , we introduced a formula for the 
growth of a Bose-Einstein condensate, in which growth 
resulted exclusively from stimulated collisions of atoms 
where one of the atoms is left in the condensate. This 
gave a very simple growth formula which predicts a rate 
of growth of the order of magnitude of that observed 
in current experiments ^^]. The direct stimulated ef- 
fect must eventually be very important, once a signifi- 
cant amount of condensate has formed, but in the initial 
stages there will of course also be a significant number of 
transitions to excited states of the condensate (quasipar- 
ticles), whose populations will then also grow. As well as 
this, there will be interactions between the condensate, 
the quasiparticles and the atomic vapor from which the 
condensate forms. This paper extends the description of 
the condensate growth to include these factors. 

All of these effects are encompassed by the descrip- 
tion given in 0j5|. However, the practical extension of 
this description to take account of the additional effects 
would involve a calculation of all the cigenfunctions for 
the trapped condensate, and the detailed summation over 
all processes involving these. Since the number of levels 
involved is of the order of tens of thousands, this could 
be a formidable task. However some quite reasonable 
estimates can be made for the overall effects of these 
processes, and from these we can derive a set of easily 
solvable differential equations for the populations of the 
condensate and the lower energy quasiparticles. 

As in our previous work, we divide the states in the po- 
tential into the condensate hand, Rc, which consists of 
the energy levels significantly affected by the presence of 
a condensate in the ground state, and the non-condensate 
band, Rnc, which contains all the remaining energy lev- 
els above the condensate band. The division between the 
two bands is taken to be at the value, E^. 

The picture we shall use assumes that Rnc consists 
of a large "bath" of atomic vapor, whose distribution 
function is given, for the energy levels greater than a 
value which we shall call E max (with .Emax > En), by a 
time-independent equilibrium Bose-Einstein distribution 



{exp[(£ - n)/k B T] - I}" 1 . The value of E max will be 
assumed to be small enough for the majority of atoms to 
have energies higher than E max , so that this part of the 
bath can be treated as being essentially undepleted by 
the process of condensate growth. 




Fig.l: Modification of the equilibrium level occupation by 
condensate growth, for a 3D harmonic trap. The cumulative 
occupation of states N(E) with energy below E, is shown for 
a gas of noninteracting atoms as a dotted line, and after con- 
densate growth (when the ground level energy has increased 
by (i(N)) is shown as a dashed line. The solid line gives 



the occupation f(E) — 



dN(E) 



after condensate growth. The 



values of fi,fj,(N), and ksT correspond to the MIT sodium 
condensate of 10 6 atoms at 1.2 fiK. 

The energy levels between Er and E m&K are taken to 
have a time dependent population, since the continuation 
of the equilibrium Bose-Einstein formula to lower ener- 
gies eventually leads to an unrealistic situation in which 
the populations and transition rates become too large 
for the populations to be considered to be constant — 
this demonstrates that the initial condition in which the 
vapor is at a positive chemical potential cannot apply 
for all energies. The choice of E max is thus determined 
as a lower limit to the equilibrium distribution, with the 
distribution in the range between E max and Er treated 
as time dependent, and computed as part of the growth 
kinetics. 

The value Er, above which the energy levels are 
taken as unperturbed, was fixed at 2/i. This value and 
the ground state energy level — the chemical potential 
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/i(iV) — put bounds on the energy levels of the states in 
between. As a simple expression of this fact, the levels 
between n(N) and Er are determined by interpolating 
linearly between the two extremes, using a density of 
states Af[E — /i(/V)] 2 , where Af is a normalization chosen 
so that the cumulative number of states matches the cor- 
responding cumulative value for an unperturbed three 
dimensional harmonic oscillator when E = Er. As an 
illustration of the effect of this, we show in Fig.l the cu- 
mulative occupation N(E), and the occupation per unit 
energy interval f(E) = dN(E)/dE, when the system is 
in equilibrium. (The condensate population itself is not 
shown.) 
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Bath of higher energy atoms 







Fig.2: The transitions being considered: Left — scattering; 
Right — Condensate growth. 

The dynamics we will consider will arise from two kinds 
of process as illustrated in Fig.2. 

Growth: A collision between a pair of atoms initially in 
the bath of atomic vapor results in one of the atoms hav- 
ing a final energy less than E max . 

Scattering: A collision between an atom initially in an 
energy level below E m&x and a bath atom transfers the 
first atom to another energy level below £ max . 

Our treatment therefore omits any scattering between 
atoms which both have energies less than £ max , which is 
reasonable if the number of atoms in the bath is almost 
100% of the total number of atoms. 

The mechanism of condensate growth, as in [jl],^) , un- 
der appropriate approximations, gives rise to equations of 
motion for the number of atoms in the condensate band 
iV, and the number of quasiparticlc excitations n m (with 
energies e m < E max ) in the condensate, which can be 
written as follows. First define 

h+ = 2W++(N) [(1 - e^-» +e ^/ k * T )n m + l] , (1) 
h- = 2W m +(N) \(1 - e(-^ N ^ +e ™V kBT )n m + l] . (2) 



Then the multilevel growth equations are 



N = 2W + 



(N) ((] 



MN)-^)/k B T 



(3) 



Here 2(N + 1)W + (N) is the transition rate for an atom 
to enter the condensate level as a result of a collision 
between two atoms in the "bath" of atomic vapor — the 
factor N + 1, takes account of both the "spontaneous" 
term, and the "stimulated" term induced by the pres- 
ence of the condensate. The reversed process occurs with 
the overall rate 2NW+{N)eS^-^/ kBT — that is with 
no "spontaneous" term, and with a factor dependent on 
the difference of the chemical potential fi of the vapor, 
and that fJ.(N) of the condensate. As a result, equilib- 
rium in the large TV limit occurs at equality of the two 
chemical potentials. 

During the process of BEC formation, the spectrum 
of eigenvalues makes a transition from the unperturbed 
spectrum of trap levels to the case where the spectrum is 
strongly affected by the condensate in the ground state. 
The Bogoliubov spectrum of a condensed gas is valid in 
the case where the number of particles in the conden- 
sate, no, is so large that it is valid to write no ~ N. 
Thus, during the initial stages of condensate formation, 
where this is not true, one must use another formalism. 
In this paper we will consider the situation in which the 
interaction between the particles is very weak, as is in 
practice the case. This means that we will be able to use 
the unperturbed spectrum for the initial stages of con- 
densation, and only use the Bogoliubov description once 
enough condensate has formed to make the effective in- 
teraction rather stronger. 

The basic formalism of || can still be carried out in 
this case, and the modification that is found is rather 
minor — essentially, we make the substitution N — > no in 
the chemical potential and the W + (N), W ++ (N) func- 
tions, and set — » 0, since this term comes from 
the mixing of creation and annihilation operators which 
arises from the Bogoliubov method. In order to simplify 
the equations we also group the levels in narrow bands 
of energy with gj- levels per group, and for simplicity use 
the same notation nu now for the number of particles in 
the energy band with mean energy e^. (This corresponds 
to the ergodic assumption used in Q.) We then deduce 
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«0 | growth = 2W + (n Q ) 
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The growth equations can be modified to include the 
terms derived in Q corresponding to the scattering of 
particles in the condensate band by the vapor particles. 
This is equivalent to scattering of particles by a heat 
bath, which leads to a rate equation for scattering of the 
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form (Where N km = 1/ (cxp[(e fc - e m )/k B T] - 1), and by 
k > m we mean e k > e m ) 



|scatt 



2_,7m/s {^VmfcnA;(n m + 5m) - (N m k + l)(n& + 5fc)™m} 
y^Jkm {(N km + l)n k (n m + g m ) - N km {n k + g k )n m } . 

(7) 



k>m 



The formulae of QKIII give precise methods for comput- 
ing the coefficients 7fc m , but we can simplify their com- 
putation by adapting the kinetic equation of Holland el 
al. ||. This methodology is based on a model in which 
the trap levels are all treated as being unaffected by the 
presence of the condensate, which should suffice as a first 
approximation. To apply it to this situation, we assume 
all the levels with energies greater than -E max are ther- 
malized, and sum out over these levels. The working is 
essentially straightforward, and yields an equation for the 
n rn variables in the form (with M the mass of the atom 
and a the scattering length) 
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where T(T) = >e e ~ e m/k B T j ias a vame which 
depends on the spectrum of energies. For an isotropic 
3-dimcnsional harmonic oscillator with frequency u>, the 
energy levels above the zero point are e„ = nhuj, so that 
we find 



r(T) 



I _ e -huj/k B T ' 



(9) 



This corresponds to essentially to (|7|) when one makes 
the correspondences 



N kn 

1km 
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1 + N k , 



1 



(10) 



7mfe 



9k 



with k > m. (11) 



The equation for both growth and scattering is now 
given by adding (||) to ( |io| ) 



growth ~t~ \ scatt i 



(12) 



where, for m = 0, we use (§ instead of (|). 

The overall evolution of the system can now be found 
from the numerical solutions to (U3). The parameters 



used were chosen so as to be in approximate agreement 
with the experimental work being conducted at MIT, 
where the growth of Bose-Einstein condensates of 23 Na is 
being studied. In contrast to the estimate in |l| in which 
the bath distribution was approximated by a Maxwell- 
Boltzmann distribution, in this computation we use the 
full Bosc-Einstcin distribution, truncated at a lower en- 
ergy of £ m ax, since lower energies are described by the 
n m variables. 

In applying the theory two major approximations are 
made. Firstly, the W+ + (N) functions were assumed 
to be equal to the W + {N) function, since the actual 
M / + + (iV) terms are not easily calculated. The justifi- 
cation for this is that the Wj+ + (7V) terms represent an 
averaging over all the levels contained in the mth group, 
and as such they may be expected to be of the same or- 
der of magnitude as W + (N) . As a validity check, it was 
found that the effect on the condensate growth rate was 
small when the W+ + (N) terms were altered by a factor 
in the range 0.5 — 2. 

In current BEC experiments the confining harmonic 
potential is normally anisotropic, whereas |6| was re- 
stricted to an isotropic trap. The second approximation 
is therefore that the scattering rate factor T(T) is equal 
to that for an isotropic 3D harmonic oscillator with fre- 
quency equal to the geometric mean frequency of the 
anisotropic trap. The precise value of this factor was 
found to have little effect on the solutions, so long as it 
was greater than about one tenth of the value given by 

The condensate rate of growth depends on the number 
of groups of levels considered in the model — modeling 
more groups in the condensate band increases the rate of 
growth, which approaches a limiting value. The number 
of groups modeled was therefore chosen as large as pos- 
sible, but it was required that there were at least 4 levels 
in the first group of levels above the condensate level. 
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Fig.3: Condensate growth for Sodium: Dotted line — Growth 
of the condensate for the uncorrected model of Q; Dashed 
line — scattering is neglected, but W + is given by using the 
Bose-Einstein correction; Solid line — with scattering and the 
Bose-Einstein correction. In all cases the initial amount of 
condensate at t = was 500 atoms. 

The initial populations for the groups in the condensate 
band were generated by a method which models the ex- 
perimental procedure. We start at some initial i « 
with the bath at a chemical potential /i w 0, and evolve 
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the equations of motion until equilibrium is reached; at 
t = we then change /i to a positive value such that 
fi = ju(iVfi na i), where iVfi na i is the final number of atoms 
in the condensate. Changing the populations at t = 
merely moves the growth curve slightly forward or back- 
ward in time, with no effect on its basic shape. We 
present a sample of the results obtained in Fig. 3. The 



trap parameters, v x = 18.5Hz, v y 



250Hz, tem- 



peratures and final condensate number are chosen in the 
range presently being investigated for sodium. 

The typical behavior of the noncondensate levels is 
shown in Fig. 4, for the case where the initial occupa- 
tion of all groups below En is chosen to be zero, in or- 
der to show the speed of the relaxation process. The 
initial growth is in the population of the noncondensate 
levels — their occupations can become very large, but this 
is because the numbers of levels in each group are very 
large, so that there is little degeneracy, i.e., the number 
of atoms per level does not significantly exceed one. Even 
in the lowest group the degeneracy is no more that about 
100. The moment the condensate achieves a significant 
degeneracy, the stimulated process takes over, and imme- 
diately draws the excess population of the noncondensate 
groups into the condensate, well before the full conden- 
sate occupation is achieved. Thus, apart from the initial 
transient, the condensate growth occurs essentially by 
the same mechanism as in [El , with the modification that 
the distribution over the noncondensate levels changes 
slowly in response to the change of the condensate chem- 
ical potential. 

If scattering is entirely neglected the populations of the 
lowest noncondensate levels become several times larger 
than those of the the condensate before settling to their 
very much lower equilibrium values. However, the inclu- 
sion of even as little as 0.1% of the strength of scattering 
used here eliminates that effect almost entirely. 
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Fig.4: Growth of noncondensate levels — the condensate num- 
ber is the almost vertical black line. 

In contrast to our work, in which explicit use is made of 
trap eigenfunctions, other descriptions of conden- 

sate growth either do not treat the case of a trapping po- 
tential, or consider only the case where trapped situation 
is not essentially different from the untrapped situation. 
Nevertheless, the equations we use have a strong con- 
nection with those of Stoof H, but their interpretation 



is different, since Stoof's Fokkcr-Planck equation treats 
the untrapped case, where the low-lying levels are ex- 
tremely closely spaced, and relative coherences between 
these levels become important. The current experiments 
are probably closer to the situation of our model, since 
the growth rate is somewhat slower than the lowest trap 
frequency. 

We can conclude from the results of this model that 
the main effect of the inclusion of the scattering and the 
computation of W + using the full Bose-Einstein formula 
is to speed up the condensate growth by up to one or- 
der of magnitude, the precise speedup depending on the 
condensate size and temperature. Even though we have 
only estimated the scattering and transition rates for the 
quasiparticle levels, the results are reasonably predictive, 
since it would be hard to credit the scattering or the W ++ 
factors as being very different from the values assumed, 
and also because the numerical predictions are not very 
sensitive to these precise values. Precise predictions will 
involve the detailed computation of these effects, rather 
than their estimation. 
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